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O Abstract. We review briefly direct and indirect methods of measuring the masses of black holes 
in galactic nuclei, and then focus attention on supermassive black holes in active nuclei, with 
special attention to results from reverberation mapping and their limitations. We find that the 
O |i intrinsic scatter in the relationship between the AGN luminosity and the broad-line region size is 
very small, ~0.11 dex, comparable to the uncertainties in the better reverberation measurements. 
We also find that the relationship between reverberation-based black hole masses and host-galaxy 
bulge luminosities also seems to have surprisingly little intrinsic scatter, ~ 0.17 dex. We note, 
however, that there are still potential systematics that could affect the overall mass calibration 
at the level of a factor of a few. 
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1. Introduction 



With the advent of suitable technology for high angular resolution spectroscopy, from 
space with Hubble Space Telescope and from the ground with adaptive optics assistance, 
it has become possible to measure the masses of the central black holes in many nearby 
t— I . galaxies. Observations of nearby galaxies have led to the identification of scaling relation- 
ships that then allow us to estimate masses of black holes in distant galaxies and thus 
•'■h . determine, at least in principle, the mass function for supermassive black holes, both 
locally and over cosmic time. Consequently, there have been tremendous advances in our 
understanding of the evolution of the supermassive black hole population and that of the 
galaxies that host them over the history of the universe. Despite this incredible progress, 
it is important to understand that supermassive black hole measurement is a field in its 
infancy: even the methods deemed most reliable have systematic uncertainties that have 
not been mitigated in a completely satisfactory way. Masses based on stellar dynamics 
and gas dynamics are highly model dependent and are only as reliable as the underly- 
ing assumptions. An obvious example is the recent work of Gebhardt & Thomas (2009) 
who found that the mass of the black hole in M 87, widely regarded as one of the most 
secure measurements, changed by a factor of two with the introduction of a dark matter 
halo into the models. It is also worth noting that only two black hole measurements, 
the Milky Way (based on orbital motion of stars in the Galactic Center) and NGC 4258 
(based on motions of megamasers) , meet the formal criterion of measuring the mass on 
a sufficiently compact scale that the black-hole nature of the central source is proven. In 
all other present measurements, it remains an article of faith that the central source is 
a singularity, although it is hard to see practically how these masses could be anything 
else. We will therefore adhere to the convention of referring to these objects as black 
holes. 

Fundamentally, mass measurements are based on how the central black holes accelerate 

1 



2 



B.M. Peterson 



nearby matter, either stars or gas. The advantage of using stellar dynamics is that stars 
respond only to gravitational forces. The corresponding disadvantage is that high angu- 
lar resolution is necessary; measurement of stellar motions at large distances from the 
nucleus makes the mass determinations more uncertain. Gas motions, on the other hand, 
allow us to probe closer to the nucleus; indeed, in the case of reverberation mapping, 
angular resolution is irrelevant since time resolution substitutes for it. The disadvantage 
of using gas motions is, of course, that gas also responds to forces other than gravity, 
with radiation pressure being the main concern (see |j5j) . 

It is important to distinguish between direct and indirect methods of mass measure- 
ment. Direct methods are based on dynamics of gas or stars accelerated by the central 
black hole. This would include stellar dynamics, gas dynamics, and reverberation map- 
ping. Indirect methods, on the other hand, are based on observables correlated with the 
mass of the central black hole. Indirect methods include the Mbr— c* and -Mbh - -^buigc 
relationships, the fundamental plane, as well as AGN scaling relationships such as the 
BLR radius-luminosity relationship that we will discuss later. 

We also sometimes refer to "primary," "secondary," and "tertiary" methods, with 
the difference based on model-dependence and assumptions required. Primary methods 
require fewer assumptions and have little model dependence: the masses of the black holes 
in the Milky Way and NGC 4258, based on proper motions and radial velocities, are thus 
certainly primary. Stellar dynamics and gas dynamics are sometimes also regarded as 
primary methods as they do not hinge on other results. On the other hand, reverberation 
mapping, as it has been practiced to date, currently depends on other "primary direct" 
methods for a zero point for the mass scale, so it is technically a "secondary method," 
though it is a still a "direct method." 

Figure 1 illustrates in flowchart format how masses are determined for both quies- 
cent and active galaxies. It is slightly oversimplified in that stellar and gas dynamics 
have been used to measure central masses of a very limited number of AGNs, as dis- 
cussed below. One-dimensional (1-d) reverberation-mapping refers to measurement of 
mean emission-line response times, i.e., reverberation mapping as it has been practiced 
to date. As noted above, this requires some external calibration of the mass scale and this 
is currently provided by assuming that the Mbh - "* relationship is the same in quiescent 
and active galaxies. With two-dimensional (2-d) reverberation mapping, we aim to pro- 
duce a velocity-delay map which should give us sufficient information to model the BLR 
dynamics and determine the central masses; model-dependence will still be something of 
an issue, as it is with stellar and gas dynamics. But this will eventually put reverberation 
mapping on a par with stellar and gas dynamics as a "primary direct" method and it will 
be possible to compare directly the Mbh-ct* relationships in quiescent and active galaxies 
rather than assume that they are the same. This may occur in the very near future: the 
first reliable detections of velocity-dependent lags were reported at this conference (see 
the papers by Bentz and Denney, these proceedings). 

In Table 1, we list published masses of a few relatively nearby AGNs for which black 
hole measurements have been made using more than one method. The most accurate mea- 
surement is provided by megamaser motions in NGC 4258. However, megamaser sources 
like this are exceedingly rare. Moreover, megamasers occur in type 2 AGNs where our 
view of the central engine is obscured, thus making reverberation mapping impossible. 
Mass measurements by stellar dynamics and/or gas dynamics have been made for a 
handful of nearby AGNs (e.g., Hicks et al., these proceedings), but in both cases, signif- 
icant future progress depends on dramatic improvements in angular resolution. Thus, if 
only by default, reverberation mapping is the near-term future for AGN black hole mass 
measurement. The downside to this is that reverberation mapping is a resource-intensive 
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Megamasers 


38.2±0.1 [1] 


N/A 


N/A 


Stellar dynamics 


33 ± 2 [2] 


7-20 [3] 


^ 70W 


Gas dynamics 


25-260 [5] 


20 +io [6] 




Reverberation 


N/A 


7.63li; 6 ^ 171 


46 ± 5 [8] 


Indirect methods: 




M BH -ai 9] 


13 


25 


6.1 


R-L scaling 1101 


N/A 


15 


29-120 



[11 Herrnstein et al. (2005). [2] Siopis et al. (2009). [3] Davies et al. (2006). [4] Onken et al. (2007). 
[5l Pastorini et al. (2007). [6] Hicks & Malkan (2008). [7] Denney et al. (2010). [8] Bentz et al. 
(2006b). [91 Giiltekin et al. (2008). [10l McGill et al. (2008). 

technique based on high-precision spectrophotometry. However, as we see below, rever- 
beration results already provide us with a shortcut to mass estimation. 



2. Reverberation-Based Mass Measurements 

Reverberation mapping (Blandford & McKec 1982; Peterson 1993) makes use of the 
time-delayed response of the broad emission lines to continuum variations to provide a 
size estimate for the broad- line region (see Peterson 2001 for a fairly thorough tutorial). 
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Figure 1. Methods and applications of black hole mass measurements in galactic nuclei. 
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At the present time, emission-line time delays, or lags r, have been measured for ~ 45 
AGNs, in nearly all cases for H/3, but in many cases for the other Balmer lines and in a 
few cases for multiple lines extending into the ultraviolet. We can then obtain the mass 
of the black hole (or, more accurately, the mass enclosed out to the distance R = cr) by 
combining this with the emission-line width AV, i.e., 

M BH = y G ' , (2.1) 

where / is a dimcnsionlcss factor of order unity that depends on the structure, dynamics, 
and orientation of the BLR. If, for example, the BLR is a thin ring of material in a 
Keplerian orbit of inclination i around the black hole, then / oc 1/ sin 2 i. Of course, the 
real geometry of the BLR must be much more complex than a simple ring or disk; unified 
models suggest that we see type 1 AGNs at inclinations 0° < i < 45°, in which case the 
observed line-of-sight velocities in a ring or disk are a small projection of the orbital 
speed. Indeed, for i < 20°, the projected line-of-sight velocity width of the line is so 
small that we would underestimate the mass of the black hole by more than an order of 
magnitude. We note in passing that emission-line lags are unaffected by inclination as long 
as the system has axial symmetry and the line emission is isotropic. Nevertheless, there is 
evidence, mostly from radio-loud systems, that inclination is an important element in the 
determination of /. For example, the ratio of core-to-lobe power in double radio sources 
correlates inversely with line width (Wills & Browne 1986); core-dominated systems, 
where we are looking down the system axis, have narrower lines than systems at larger 
inclination. Similarly, flat-spectrum radio sources, in which again we are observing close 
to the axis of the system, have narrower emission lines than steep spectrum sources, 
which are seen at higher inclination (Jarvis & McLure 2006). But the differences are 
not extraordinarily large: for example, Jarvis & McLure find that for their sample of 
radio-loud AGNs, sources with a > 0.5 have a mean FWHM of 6464 km s -1 for H/3 and 
Mgn, while those with a < 0.5 have a mean width of 4990 km s -1 . These results clearly 
indicate that (1) inclination is important, but (2) the BLR gas has considerable velocity 
dispersion in the polar direction. Given this, we now have to recognize that a direct 
comparison between reverberation-based masses and, say, stellar or gas dynamics, will 
depend very much on the generally unknown inclination of the BLR. We should generally 
not expect good agreement between reverberation masses and masses measured by other 
techniques for individual galaxies unless the inclination is known or strongly constrained. 

On the other hand, given a sample of AGNs, we can determine an average value (/) by 
normalizing the AGN Mbh-ct* relationship to that of quiescent galaxies. A first attempt 
to do this yielded a value (/) = 5.5 ± 1.8 (Onken et al. 2004): while this intuitively 
seems a little high, we can again plausibly attribute this to inclination effects and our 
predisposition the inclinations of type 1 AGNs are typically rather small. The main 
liability of this particular calibration of the mass scale is that it is based on low-redshift, 
low-luminosity objects. All of the a* measurements are based on the Can triplet which 
is unfortunately redshifted into atmospheric water vapor bands at z > 0.06. We are now 
attempting to measure at higher rcdshift by using observations of the CO bandhead in 
the //"-band (1.6 /zm), a process greatly aided by use of adaptive optics with integral field 
units that concentrate the bright AGN into the central few pixels and at the same time 
integrate a few arcseconds in each coordinate without loss of spectral resolution (see the 
papers by Dasyra and Grier in these proceedings). In addition to new results at higher 
redshifts and luminosities, the LAMP collaboration (Bentz, these proceedings) is adding 
additional masses and velocity dispersions at the low-mass end of the AGN distribution. 
As a consequence of expanding the range in black hole mass, the slope of the AGN Mbh~ 
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<7* relationship is becoming better denned and shows improved consistency with that for 
quiescent galaxies (Woo, these proceedings). 

If we assume that the AGN Mbh-c* relationship has the same zero point and slope 
as that for quiescent galaxies, a maximum likelihood analysis places an upper limit on 
the intrinsic scatter in the relationship of AlogMeH ~ 0.40 dex, which is consistent with 
what is found for quiescent galaxies (Giiltekin, these proceedings). This is a reassuring 
indication that the masses derived from reverberation mapping are reliable. There are, 
of course, other such indications, including the direct comparisons of a few cases as in 
Tabic 1. Also, in each case in which lags have been measured for multiple emission lines 
in a single source, there is an anticorrelation between lag and line width that is consistent 
with a constant virial product AV 2 R (Peterson & Wandel 1999). Finally, AGNs show 
the same correlation between black hole mass and bulge luminosity (the MBH _ ^buige 
relationship) that is seen in quiescent galaxies (Bentz et al. 2009b). Indeed, a maximum 
likelihood analysis gives an upper limit to intrinsic scatter, AlogMeH ~ 0.17 dex, which 
is actually smaller than the intrinsic scatter in this relationship for quiescent galaxies 
(A log Mbh « 0.38 dex; Giiltekin et al. 2008). 



3. BLR Scaling with Luminosity: The R—L Relationship 

The emission-line spectrum from an ionized gas, setting aside elemental abundances, 
is largely controlled by the particle density n and an ionization parameter 

u- gro (31) 

where Q(H) is the number of hydrogen- ionizing photons emitted per second by the central 
object. To some low order of approximation (principally ignoring the Baldwin Effect), all 
AGN spectra have similar emission-line flux ratios and similar emission-line equivalent 
widths, suggesting that U and n do not vary much from one source to another. If we 
further assume that Q(H) oc L where L is the luminosity of the central source in some 
arbitrary band, we are led to the naive prediction that R oc L 1 / 2 . 

Thus since the early days of photoionization equilibrium calculations, some such re- 
lationship between the BLR radius and the continuum luminosity (hereafter simply the 
R-L relationship) of the power-law form R oc L a had been anticipated. Observationally 
the R-L relationship for H/3 first emerged as statistically significant with a slope a w 0.7 
with the addition of 17 PG QSOs by Kaspi et al. (2000) to the existing similar-sized 
reverberation-mapping database on Seyfert 1 galaxies. It is sometimes stated that this 
result was enabled by extension of the luminosity range by nearly two orders of magni- 
tude; it is more correct to say, however, that it was due to the fact that the QSOs are so 
much more luminous than their host galaxies that starlight contributes negligibly to the 
luminosity measurement. Indeed, it is surprising how much starlight contaminates AGN 
luminosity measurements in the optical, though this problem is certainly exacerbated in 
reverberation monitoring data as large spectrograph apertures and spectral extraction 
windows are used to mitigate the effects of variable seeing on the total flux measurements. 
Bentz et al. (2006a, 2009a) have attempted to account for starlight contamination of op- 
tical luminosities of AGNs by using unsaturated high-resolution optical images obtained 
with HST/ ACS. The surface brightness distribution is modeled so that the AGN point 
source can be removed and the galaxy contamination estimated by simulated aperture 
photometry When this is done, the slope of the R-L relationship for H/3 is found to be 
a w 0.5, remarkably consistent with the naive prediction. 

This leads us to the remarkable realization that it is possible to estimate masses of 
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the black holes in AGNs from only a single spectrum: by measuring the luminosity, we 
infer the BLR radius and we can combine this with the emission-line width to measure 
the mass (Wandel, Peterson, & Malkan 1999). This is variously referred to as the "pho- 
toionization method," "single-epoch spectrum method," or simply "masses from scaling 
relationships." 

An interesting question at this point is how much intrinsic scatter is in the R-L 
relationship? This is of keen interest for two reasons: first, this will set a fundamental limit 
on the accuracy with which we can estimate masses from single spectrum, and second, it 
dictates future observing strategies for further refinement of the R-L relationship. If the 
intrinsic scatter is large, then many more reverberation measurements will be needed to 
beat down the statistical noise in the relationship. On the other hand, if the scatter is 
small, then improvements in the determination of the R-L relationship will come from 
obtaining better reverberation data on a smaller number of objects. 

To proceed with this, we want to start by using only the very best data sets, since 
these are ones where intrinsic scatter is most easily isolated: there remain elements of 
both art and luck in reverberation studies, and the ugly truth is that not all of the light 
curves give as clean a result as we would like. To minimize the impact of the lower-quality 
data, Catherine Grier and I independently visually inspected all available light curves 
(for the optical continuum and H/3 emission line only) with the intent of identifying the 
"best" reverberation data: we used a very simple criterion, that you could easily see the 
same pattern of variability in the continuum and H/3 light curves and could estimate 
the lag simply by inspection (see Figure [2]). About half of the existing light curves met 
this criterion. A maximum likelihood analysis of these data indicates that the intrinsic 
scatter amounts to A logi? « 0.11 dex. This is really quite remarkable when one considers 
that the typical formal errors on this subset of the best reverberation data amount to 
A log R « 0.09 dex. From this, we conclude that, at least over the luminosity and redshift 
range where it has been calibrated (41.5 < logALA(erg s _1 ) < 45 at A = 5100A and 
z w 0), the R-L relationship is as statistically effective as reverberation for obtaining 
BLR sizes and central black hole masses. If you are wondering why we emphasize the 
"statistical" aspect of R-L scaling, we refer you to the results of the "indirect methods" 
in Table 1: the results for individual sources can be quite misleading, as they certainly 
are in the case of NGC 4151, which is a notorious outlier in the Mbh _ c* relationship. 



4. The R—L Relationship at High Redshift and Indirect Mass 
Measurements 

4.1. The R-L Relationship for UV Lines 

With the R-L relationship, we are able to explore the black hole mass function, not 
only locally but at high redshift, enabling us to trace the history of black hole growth. 
Some exploratory work has been done on this and in fact there are claims that the 
Mbh~o~* and MeH-^buige relationships evolve over time (e.g., Woo, these proceedings), 
although at least some published claims of evolution of the Mbh _ o* relationship are 
clearly attributable to Malmquist bias. A lot of additional careful work will be required 
to sort this out. 

An obvious problem with estimating black hole masses at high redshift, of course, is 
that the H/3 emission line is redshifted out of the visible window at only modest redshift. 
We are thus forced to use other emission lines for which reverberation measurements are 
actually quite scarce: of the other potential broad lines for mass measurement, there is 
but a single reliable measurement of a Mgn lag (Metzroth, Onken, & Peterson 2006) 
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and a bare handful of C iv lags, though these span quite a large range in luminosity, 
39.5 < logAL A (erg s" 1 ) < 47 at A = 1350A in the rest frame (Kaspi et al. 2007 and 
references therein). One of these objects is at z = 2.17, and the rest are at z < 0.06. 
The Civ A1549 emission line is especially important as it allows us to reach quite large 



rcdshifts in the optical. Vestergaard (2002) first used Civ to show that quasars with 
black hole masses of ~ 10 a M Q were already assembled by z ~ 4. She assumed that the 
C IV R—L relationship has the same slope as H/3 and calibrated the C iv-based mass scale 
by using single-epoch spectra of reverberation-mapped AGNs. 

4.2. Characterizing the Velocity Field 

If there is an "Achilles' heel" in determining masses based on individual spectra, it is 
probably not the R-L relationship, but rather how we characterize the velocity field of 
the BLR. The line profile is expected to be sensitive to both inclination and Eddington 
rate (e.g., Collin et al. 2006). 

In reverberation mapping experiments, we obtain the best results by measuring the 
line width in the "rms spectrum" formed by computing the pixel-by-pixel variance in 
all of the spectra obtained in the reverberation experiment. This isolates the variable 
part of the emission line, which arises in the very gas for which we are measuring the 
time delay. There are two common ways to characterize the line width, either FWHM or 
the line dispersion erg, which is the second moment of the line. The latter seems to give 
somewhat better results (Collin et al. 2006; Peterson et al. 2004). For either line- width 
measure, the modest scatter, e.g., in the -MBH - £buigc relationship, suggests that the mass 
measurements are not particularly sensitive to how we characterize the line width and 
thus the velocity field of the BLR. 
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Figure 2. The R-L relationship for H/3. The luminosity is ALa(5100 A) and the BLR radius is 
measured in light days. The open circles indicate the highest-quality measurements. The slope of 
the fit to the highest- quality data is indistinguishable from that to all the data, but the scatter 
is only 0.11 dex. 
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The situation is rather different with singie-epoch spectra where contamination by 



other features becomes a serious issue. Denney et al. (2009) have explored this in some 



detail for H/3 by seeing how well single spectra from reverberation campaigns reproduce 
the actual reverberation-based masses. They find some obvious effects (e.g., sensitivity of 
FWHM to proper removal of the narrow-line component) and some more subtle effects 
(e.g., a i is badly affected by blending, particularly when the AGN continuum is in a faint 
state and the lines are especially broad). 

Arguments about whether particular UV lines can be used to estimate masses are 
really about how well the widths of different lines in AGN spectra are correlated. In the 
case of Mgn A2798, the width of the Mgn line appears to be a good surrogate for the 
width of H/3 (Woo, these proceedings). Shen et al. (2008), using a large collection of 
SDSS spectra, find that log ([FWHM (H/3)] / [FWHM(Mgn)]) = 0.0062, with scatter of 
only ~ 0.11 dex. Onken & Kollmeier (2008) find that the masses derived from H/3 and 
Mgn differ as a function of Eddington ratio, although they argue that it is possible to 
correct for this bias. 

The accuracy that can be obtained using the C IV line remains somewhat controver- 
sial (see the contributions to these proceedings by Vestergaard, Netzer, and Trakhten- 
brot) . On the positive side, the limited existing reverberation data for C iv suggest that 
the R-L relationship for Civ has the same slope as H/3 (Kaspi et al. 2007). Moreover, 
reverberation-based masses based on C IV are consistent with those based on every other 
line (i.e., the virial product AV 2 R/G is consistent for all the lines). On the other hand, 
since C iv is a resonance line, we often see absorption in its blue wing due to outflows, 
complicating accurate measurement of the line width. Indeed, we sometimes see, par- 
ticularly in the case of narrow-line Seyfert Is (NLSls), an extended blue wing of the 
emission line, again presumably due to outflowing gas, which makes measuring the line 
width problematic. It seems likely, however, that we should be able to calibrate out 
luminosity-dependent effects and simply avoid spectra with strong absorption features 
such as BALs. 

4.3. Use of Scaling Relationships 

The scaling relationships that lead to black-hole mass estimates must not be used blindly; 
indeed they should be used with great caution. We need to keep in mind that when we 
think we are measuring mass, we are really measuring 

M BH oc AV 2 R oc AV 2 L 1/2 . (4.1) 

Similarly, when we think we are measuring Eddington ratio, we are really measuring 

L L L L 1 / 2 

L Edd a M B „ K A^ K AF' (4 ' 2) 

It is important to keep in mind that any correlations among mass, Eddington ratio, and 
luminosity must be eyed with great suspicion. 



5. On the Possible Importance of Radiation Pressure 

Given a nominal bolometric correction, most AGNs are thought to be radiating at 
about O.lLEdd where ^Edd is the Eddington luminosity, the maximum luminosity at 
which the source is stable against disruption by radiation pressure. In the case of objects 
like NLSls, the luminosities may approach L^dd- It stands to reason that radiation 
pressure on the BLR gas might be an important dynamical force. In an elegant paper, 
Marconi et al. (2008) argue that radiation pressure could be an important factor in the 



Central Black Hole Masses 



9 



BLR dynamics. Radiation pressure, like gravity, is diluted geometrically (i.e., oc R~ 2 ), 
making it hard to distinguish between a high central mass plus radiation pressure and 
a lower central mass and no radiation pressure: failure to account for radiation pressure 
could thus lead us to underestimate the mass of the black hole. Marconi et al. argue that 
the relative importance of radiation pressure ultimately comes down to the inertia of the 
BLR gas clouds: if the clouds have sufficiently high column density (iVn > 10 23 cm~ 2 ), 
then the effect of radiation pressure on the gas dynamics is negligible. Interestingly, the 
most successful photoionization equilibrium calculations suggest column densities of this 
order. Marconi et al. then modify equation (|2.1[) to a form like 

M BH = fAV 2 R/G + gL, (5.1) 

where the first term on the right hand side is the mass based on the R-L scaling rela- 
tionship and the second term accounts for radiation pressure. Marconi et al. claim two 
successes with this formulation: first, the scatter between single-epoch mass estimates 
and reverberation measurements is reduced and the difference is no longer luminosity- 
dependent (their Figure 1), and second, NLSls no longer fall below the normal Mbh-ct* 
relationship. It is not clear, however, that anything physical can be inferred from this. The 
additional luminosity dependence of equation (|5.1j) can be attributed to the fact that the 
reverberation-based sample of AGNs is strongly Malmquist biased: the reverberation- 
mapped AGNs were selected because they are among the apparently brightest known 
AGNs, and there is a clear correlation between their black hole masses and luminosities. 
Thus the second term in equation (|5.1|) is bound to reduce scatter, regardless of whether 
or not radiation pressure is actually physically important. Moreover, the argument about 
NLSls appears to be circular: if accounting for radiation pressure increases their black 
hole masses enough to place them on the Mbh _ ct* relationship, then their Eddington 
rates correspondingly drop down to those comparable to other AGNs, which then begs 
the question of the origin of the unusual properties of NLSls. Also, Netzer (2009) drew 
large samples of Type 1 and Type 2 AGNs from SDSS and showed that while their [O in] 
luminosity distributions and conventionally computed black hole masses have similar dis- 
tributions, the radiation-pressure corrected mass distributions were markedly different, 
arguing against the importance of the radiation pressure term in equation (I5.ip . Marconi 
et al. (2009) countered that Netzer's result is an evitable consequence of assuming a 
fixed column density for the BLR clouds; unfortunately, this argument exposes the fact 
that there is currently no accurate formulation for estimating black hole masses in the 
presence of radiation pressure. 

At the present time, whether or not radiation pressure needs to be accounted for in 
our mass calculations is not clear, though many of us continue to work on this issue. We 
believe that NLSls probably provide the best testing ground, and one obvious step is to 
further investigate the Mbh _ c* relationship in these sources. Much of the previous work 
on NLSls has necessitated using the [O in] narrow line widths as a surrogate for stellar 
velocity dispersions, and this is one area where we can make an improvement fairly easily 
by measuring stellar velocity dispersions in the near IR. 

6. Conclusions 

While great progress has been made in measuring the masses of the central black holes 
in both active and quiescent galactic nuclei, there are still significant uncertainties. My 
own guess is that masses measured by all direct methods are probably uncertain by a 
factor of a few, though this could be an underestimate for reverberation-based masses if, 
for example, radiation pressure turns out to be important. Aside from this possibility, the 
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most significant unknown in reverberation-based masses is the inclination of the system. 
This, of course, is merely one aspect of a rather larger issue: it is also not clear that we 
have identified an unbiased way to characterize the emission-line widths that enter into 
the mass determination. The size of the BLR, at least as measured with H/3, seems to 
be remarkably well-characterized (though we could be fooling ourselves there, too) with 
uncertainties as small as A log R « 0.1 dex. 
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